HEAT KERNEL EXPANSIONS ON THE INTEGERS 



F. ALBERTO GRUNBAUM AND PLAMEN ILIEV 

Abstract. In the case of the heat equation ut = u xx + Vu on the real 
line there are some remarkable potentials V for which the asymptotic 
expansion of the fundamental solution becomes a finite sum and gives 
an exact formula. 

We show that a similar phenomenon holds when one replaces the 
real line by the integers. In this case the second derivative is replaced 
by the second difference operator Lo- We show if L denotes the re- 
sult of applying a finite number of Darboux transformations to Lo then 
the fundamental solution of ut = Lu is given by a finite sum of terms 
involving the Bessel function / of imaginary argument. 



1. Heat kernel expansions 

The subject of heat kernel expansions on Riemannian manifolds with or 
without boundaries cuts across a number of branches of mathematics and 
serves as an interesting playground for a whole array of interactions with 
physics. 



It suffices to mention, for instance, the work of M. Kac [18] on the issue of 
recovering the shape of a drum from its pure tones, as well as the suggestion 
by H. P. McKean and I. Singer [2C] that one should be able to find a "heat 



equation proof" of the index theorem. For an updated account, see [||, 23 1 . 

In the simpler case of the whole real line, the fundamental solution of the 
equation 

u t = u xx + V(x)u, u(x,0) = 5 y (x) 

admits an asymptotic expansion valid for small t and x close to y, in the 
form 



u(x,y,t) 



n=l 



When V is taken to be a potential such that L = (d/dx) 2 + V belongs to a 
rank one bispectral ring, then something remarkable happens, namely this 
expansion gives rise to an exact formula consisting of a finite number of 
terms and valid for all x, y as well as all t. For a few examples of this see 
[15 1, as well as references in Q. See also [24] where the author points to 
[21 for useful connections between the heat equation and KdV invariants. 
For a general discussion of the bispectral problem see [|io|[ and \Wk and for a 
sample of problems touching upon this area see 0, H f|0, [9|, 
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In the case above the potentials T^'s are the rational solutions of the KdV 
equation decaying at infinity, see p|. They can be obtained by a finite 
number of applications of the Darboux process starting from the trivial 
potential V = 0, see [Qj. The remarks above rest on the basic fact that the 
Darboux process maps operators of the form d 2 /dx 2 + V into themselves. 
Following some preliminary explorations in |14] we would like to use this 



"soliton technology" as a tool for the discovery of the general form of a heat 
kernel expansion on the integers. In this case, there seems to be no general 
theory predicting even the existence of asymptotic expansions. 

If one replaces the real line by the integers and looks for the fundamental 
solution of 

ut(n, t) = u(n + 1, t) — 2u(n, t) + u(n — 1, t) = Lou(n, t) 

with u(n, 0) = S nm , one obtains, in terms of the Bessel function I n {t) of 
imaginary argument, the well-known result 

u(n,t) = e- 2t I n . m {2t). 



A nice reference for this is |12[ , volume 2. In the context of I (Z) it is 
simplest to study Lq by Fourier methods and obtain for the fundamental 
solution of the heat equation the expression 



-2t 



c t(x+x- 1 ) x n-m ( ^ x _ ^.1) 



2m 

One can replace the second difference operator Lq above by an appropri- 
ate perturbation of it and look at the corresponding heat equation and its 
fundamental solution. The purpose of this paper is to describe the result 
when the ordinary second difference operator is subject to a finite number 
of Darboux factorization steps. In this case the spectrum is a finite interval 



and the factorization can be performed at each end (see (2.1)). Formula 
( |1.1| ) above will be modified properly, in (fT3|), when Lq is subject to a finite 
number of applications of the Darboux process. 

We close this introduction with the simplest nontrivial example and then 
we describe the organization of the paper. 

Let us write the operator Lq in the form Lq = A — 2Id + A -1 , where A 
stands for the customary shift operator, acting on functions of a discrete 
variable n £ Z by A/(n) = f(n + 1). We can factorize Lq as 

Lq = PqQq, (1.2) 

where the operators Pq and Qq are given by 

P = Id - ^^A- 1 and Q = A - ^±^Id, (1.3) 

with T n (5) = n + 5. Applying one Darboux step with parameter 5 to the 
operator Lq amounts to producing a new operator L\ q by exchanging the 



order of the factors in (1.2), i.e. 



L ,o - QoPa = A - ( 2 + 1 J Id + ^y*> A-i. (1.4) 

T n (5)T n+1 (5) J T n (5) 2 
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The fundamental solution to 

u t = L lfi u (1.5) 

is given by 



-2t 



u(n, m, t) 



T~rnT~n+lIn—in(2t^) tl n — m (2t) tl n — m+ l(2t) . (1.6) 



Indeed, it is obvious that u(n, m, 0) = 5 nm , and equation (|1.5|) can be verified 



using well known identities satisfied by the Bessel functions (see (4.2a) and 

It is important to point out that, when the real line is replaced by the 
integers, the operators obtained from Lq by the Darboux process are no 
longer of the form Lq plus a potential. 

The operators L, obtained from Lq by a finite number of applications 
of the Darboux process at the ends of the spectrum, have been recently 



determined, see [16, 17]. These operators belong to a rank-one commutative 
ring Ay of difference operators with unicursal spectral curve. The common 
eigenfunction p n {x) to all operators of Ay, with spectral parameter x, is also 
an eigenfunction to a rank-one commutative ring of differential operators in x 
with solitonic spectral curve, i.e. we have a difference-differential bispectral 
situation. Moreover, the functions p n {x) satisfy an orthogonality relation 
on the circle. These results will be summarized in the next section and put 
to use in later ones. 

The ring of operators Ay is also discussed in section 2 and certain opera- 
tors in it are exhibited in section 3. In section 4 we collect some properties 
of Bessel functions and then all of this is used in section 5 to obtain our 



main result, Theorem 5.1. We close the paper with one example to illus- 
trate all the steps of the proof of our main result what covers all cases when 
L belongs to a rank one bispectral ring. 



2. Rank-one bispectral second-order difference operators 

Denote by A and A, respectively, the customary shift and difference op- 
erators, acting on functions of a discrete variable n £ Z by 

A/(n) = f(n + 1) and A/(n) = f(n + 1) - /(n) = (A - Id)/(n). 

The formal adjoint to an operator 

j 

is defined to be 

3 3 

In this section we describe the operator L^s obtained by successive Dar- 
boux transformations from the operator Lq = A — 2Id + A~ 1 . The symmetric 
operator Lq admits (as d 2 /dx 2 did in the case of the real line) a unique self- 
adjoint extension in Z 2 (Z) ("limit point" case of Weyl's classification at both 
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end points). Its spectrum is the interval (—4,0). The steps of the Darboux 
transformations are as follows 

Lq = PqQo rx Li ;0 = QqPo = P\Qx rx ■ ■ ■ rx L Rfi = Qr-\Pr-i 
L R>0 + 4Id = PrQr rx L RA + 4Id = QrPr = Pr+iQr+i rx ■ ■ ■ 

rx L R>S + 4Id = Qr+s-iPr+S-1- (2.1) 



At each step we perform a lower-upper factorization, as we did in ( |1.3| ), 
of the corresponding operator and then we produce a new operator by in- 
terchanging the factors. The operator L R} $ is obtained by performing R 
Darboux steps at one end of the spectrum of Lq and S steps at the other 
end. 

For details of the exposition that follows we refer the reader to 17]. 



Let e(k, A) be the linear functional acting on a function g{z) by the for- 
mula 

(e(k,\),g)=gW(\), A e C, k > 0. 
We shall denote by Exp(n; r, z) the exponential function 



Exp(n; r, z) = (1 + z) n exp ^ 



oo 

r i z i 



where r = (r*i , r 2 , ■ ■ ■ ) ■ Let us introduce also the functions 

SUn; r) = h e (J, e - 1), Exp(n; r, z)), (2.2) 

r- 

and write 

4*j{n;r) = S^^n + j - l;r) and i/jj{n;r) = S^^n + j - l;r). (2.3) 

When e = 1, the functions Sj(n;r) are a shifted version of the classical 
elementary Schur polynomials, defined by exp ( Yl'jLi r j z ^) = Yl'jLo Sj{r)z^: 

Sj(n;r) = Sj(n +n,r 2 - n/2,r 3 + n/3, ....). 

For e = —1, the functions SJ {n; r) are of the form 

/ °° . \ 

S- 1 (n; r) = (a polynomial in n, ri,r 2 , . . . ) x (—1)" exp 2) J Vj . 

Finally, let us define 

r(n;r) =Wr A (0i(n; r), . . . ,(/> R (n;r),ipi(n;r), . . . ,ip s (n;r)) 

x(-ir 5 exp^-5^r l (-2)^, (2.4) 

where Wr A denotes the discrete Wronskian with respect to the variable n: 
Wr A (/x(n), f 2 (n), . . . , f k (n)) = det(A i " 1 /,(n)) 1 < ij < fc . 
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The purpose of the exponential factor in Q2.4D is to cancel the exponential 
factor that comes out from the functions ipj(n; r), 1 < j < S. With this nor- 
malization, the function r(n; r) becomes a quasipolynomial in the variables 
n ; r i ; r 2, ■ ■ ■ , i-e. in general, r(n, r) depends on infinitely many variables, but 
there exists a positive integer N, such that r(n, r) is a polynomial in every 
variable of degree at most N. The operator 111 (2-1) can be expressed 
in terms of the function r(n, r) via the formula 

/ 9 r(n+l;r) ^ r(n - 1; r)r(n + 1; r) t 

£fl,S = A + -2 + — log — IdH -2 A . 

\ on T(n;r) J T{n;r) z 

(2.5) 

From ( |2.4D it follows immediately that r(n, r) is an adelic tau function^ of 
the AKP hierarchy defined by R+S one-point conditions, with R conditions 
at the point zero and S conditions at the point —2. The reason why the 
end points of the spectrum of Lq, and —4, are replaced by and —2 can 
be traced back to ( |2.2| ) and Q2.3| ). After a suitable linear change of time 
variables ri,r2,... the tridiagonal operator solves the Toda lattice 

hierarchy, i.e. 

^j = m + , L ], j = i,2,... 

where (£-')+ denotes the positive difference partf] of the operator IJ ', and 
{r^}^! are related to {rj}j? =1 via linear transformation of the form 

oo 

r j = r 'j + c J fcr fc' 
fc=i+l 



see [16 1 for explicit formulas. The wave function w(n;r,z) and the adjoint 



wave function w*(n; r, z) are defined by 

. r(n; r — \z~ 1 ]) _ . . 
io(n; i, <zj = p r bxp(n; r, 




Exp(n;r, z), (2-6) 



and 



*/ \ T{n;r + [z : ]) _ x 

w (n;r, z) = bxp (n;r,z) 

T(n; r) 



oo 

l + ^2w*(n;r)z- j J Exp" 1 ^; r, z), (2.7) 



x By an adelic tau function we mean a tau function built from a plane belonging to 
Wilson's adelic Grassmannian |27j via the construction in |lq ]. 

2 Equivalently, if we think of U as an infinite matrix, {L 3 ) + is the upper part of this 
matrix, including the main diagonal. 
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where [z] = (z, z 2 /2, z 3 /3, . . . ). The wave operator W(n; r) and the adjoint 
wave operator (W~ 1 )* can be written in the form 

oo 

W{n;r) = 1 + ^ Wj (n; r)A~ j = Q(A - Id) _jR (A + Id) -5 , (2.8) 
i=i 

oo 

(W _1 )*(n; r) = l + J2 w *j( n + r ) A * ~ J = P *( A * ~ Id) _R (A* + Id) -5 , 

(2.9) 

where P and Q are finite-band forward difference operators of order R + S 
satisfying 

PQ = (A - Id) 2R (A + Id) 25 , (2.10) 

see [17, Theorem 3.3]. 

The tau function r(n; r) can be viewed as an infinite sequence (indexed 
by n) of tau functions of the standard KP hierarchy, associated with a flag 
of nested subspaces 

V : • • • C V n+X C V n C K_i C • • • , (2.11) 

with 

V n = Span{w(n;0,z),w(n + 1; 0, z), w(n + 2;0,z), ...}. (2.12) 
The plane Vq belongs to Wilson's adelic Grassmannian p7j. Let x = z + 1, 



Pn(^) = r)x ra = w(n, r, x — 1) exp — rj(x — 1) J , (2.13) 

V i=i / 

and similarly 

p* n (x) = W*(n - l;r)x- n = w*(n,r,x- 1) exp ^Tjjix - l) j j . (2.14) 

From (^) it follows easily that p n (x) are eigenfunctions of the operator 
Lr,s with eigenvalue x — 2 + a; , i.e. 

L R ,sPn(x) = (x-2 + x- l ) Pn {x). (2.15) 

Moreover, p n (x) are also eigenfunctions of a differential operator in x (with 
coefficients independent of n) and thus provide a difference-differential ana- 
log of the rank-one solutions (the KdV family) of the bispectral problem 
considered by Duistermaat and Grunbaum in flofl . 

Below we discuss the spectral curve and the common eigenfunction of 
the maximal commutative ring of difference operators Ay containing L^s- 
The correspondence between commutative rings of difference operators and 
curves was studied by van Moerbeke and Mumford [26] and Krichever |H 



(see also p2{ where the case of singular curves was treated very completely) . 
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Let us introduce the ring Ay of Laurent polynomials in x that preserve 
the flag V: 

Ay = {f(x) £ C[x,x _1 ] : f(x)V n C V n+ k, for some k £ Z and Vn G Z}. 

(2.16) 

For each /(x) G Ay there exists a finite band operator Lf, such that 

Lfp n {x) = f(x)p n (x). 

Moreover, if f(x) = ^j=mi a,jX J with a mi / and a m2 ^ 0, then the 
operator Lf has support [mi, 7712], i.e. 

JT12 

j=m 1 

We can think of Lf as an operator obtained by a Darboux transformation 
from the constant coefficient operator /(A). The ring Ay is generated by 
the functions^] 

If 1\ , 1 ( X - 1)2^+1(3; + 1)25+1 

" = 2 [ X + x J and V = ¥^sTT x R+s+i > ( 2 - 17 ) 

i.e. Ay = C[w,v]. Thus the operator L^s constructed above belongs to a 
maximal rank-one commutative ring of difference operators Ay isomorphic 
to the ring of Laurent polynomials Ay. The spectral curve of Ay is 

Spec(^v) : v 2 = (w - l) 2R+1 (w + 1) 25+1 . 

The operators Q and P* in formulas ( |2.8| ) and ( |2.9D , can be expressed in 
terms of the functions {4>i(n; r)}fL± and {ipj(n; r)}* =1 as 

Qf( n ) - Wr AOi( re i r )i • • ■,^R(n;r),ip 1 (n;r), . . . ,ip s (n;r), f(n)) 
WrA(4>i(n;r),...,(j) R (n;r),tlj 1 (n;r),...,ips(n;r)) 

and 

p * fM = Wr A * (</>f (n; r), . . . , <j>* R {n\ r),ipf(n; r), . . . , ip* s (n; r), f(n)) 

(2.19) 

with <f>*(n;t) = 4>i{n + R + S;t), 1 < i < R, and ipj(n;t) = tpj (n + R + S; t) , 
1 < j < S. Using these explicit formulas for Q and P, one can check that 

Pnix- 1 ) = T{n + 1 \ r) xp* n+1 (x). (2.20) 

Finally, (see 17, Theorem 5.2]) one can prove the following orthogonality 
relation 

^- I Pn(x)p m {x- 1 )— = T( " + If^ nm, Vn, m € Z, (2.21) 
irrz j c x T(n; t) 

where C is any positively oriented simple closed contour surrounding the 
origin, not passing through the points x = dbl. Indeed, p n (x) are rational 
functions on the Riemann sphere with poles only at x = 0, ±l,oo. The 



3 Notice that these generators differ slightly from the ones chosen in |l7 



s 
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spectral curve Spec(*4v) has cusps at ±1, hence using the fact that the 
residue of a regular differential at a cusp is always zero (see |2f| ) , we deduce 
that 

Tes x= ± 1 p n (x)p* m (x)dx = 0. 



The proof of ( ^2l|) now follows from the discrete Kadomtsev-Petviashvili 



bilinear identities and the relation between the wave and the adjoint wave 
functions in ( [2.20 ). 



3. Certain Laurent polynomials in Ay 



The main result in this section is Proposition |3.1| below which allows us 
to construct some "universal" Laurent polynomials in Ay. This result is 



crucial for the proof of Theorem 5T . As we shall see in Section 5, applying 
Proposition 3^, we can reduce the infinite sum in the computation of the 
fundamental solution of (|5.1[) - (|5.2|) to a finite sum, modulo some identities 
among the Bessel functions to be discussed in the next section. 

Proposition 3.1. Let T be > m&x(R,S), and let sq, s\, . . . , st be distinct 
nonzero integers, such that Sj = st (mod 2) and Sj + Sk ^ 0, for < j, k < 
T. Then 

T 

^k 

F S0 ,..., ST ( X ) = Yl - n X (^ — e Av - (3 - 1} 



The proof of Proposition 3.1 is based on two simple lemmas related to 
the Lagrange interpolation polynomial and the Chebyshev polynomials of 
the second kind. 

Lemma 3.2. Let q(n) be an odd polynomial in n of degree at most 2T — 1, 
and let so, s\, . . . , st be distinct positive numbers. Then 

E-FT^ 2T = 0- (3-2) 



Proof of Lemma 3.L Consider the Lagrange interpolation polynomial for 
q(n) at the nodes ±so, ±si, • • • , ±st- Since degg(ra) < 2T — 1, the Lagrange 
polynomial must be identically equal to q(n), i.e. we have 

k=0 K j^k k 3 fc=0 K j^k k J 

Hence, the coefficient of n 2T+1 on left-hand side must be zero, which gives 



(» □ 

The Chebyshev polynomials of the second kind are defined by the follow- 
ing three term recurrence relation 

2wU n (w) =U n+1 H +Z7 n _i H, n = l,2,..., (3.3) 
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with Uo(w) = 1 and U\{w) = 2w. From the last formula one sees imme- 
diately that U n (w) is an even/odd polynomial if n is an even/odd positive 
integer, respectively. 

The Chebyshev polynomials of the second kind can be obtained from the 
Jacobi polynomials by taking a = (3 = ^: 

U n (w) = (n + 1) 2 F 1 ( '3 |^— 

e p^i)! ( n (n + 1)2 - j2) ) {w - i)k - (3 - 4) 

These explicit expressions for U n (w) will be needed below. 

Lemma 3.3. If so, s\, . . . , st are distinct positive integers, then the polyno- 
mial 

T 

Qs ,...,s T (w) = ^ (3.5) 

fc=0 Sfc HjVfcl s fc s jJ 

is divisible by (w — 1) T . Moreover if sq = s\ = • • • = sx (mod 2), then 

(w 2 -lf/Q S(h ..., ST (w). (3.6) 



Proof of Lemma 3.5. From (3.4) we see that 

i.e. U£li(l) is an odd polynomial in n of degree 2k + 1. From Lemma 3.2 
it follows that 

Qig... )ST (l) = 0, fc = 0,l,...,T-l, 

hence (w-l) T ' /Q S0 ,..., ST (w). If s = si = ■ ■ ■ = s T (mod 2), then Q S0 ,..., St (vj) 
is either even or odd polynomial, which proves ( |3.6[) . □ 



Proof of Proposition \3.J\ . Notice that 

F So ,...,s k _ 1 ,s k ,s k+1 ,...,s T (x) — Fs ,...,s k _ 1 ,~s k ,s k+1 ,...,s T (x) 

= — rf T-2 2T^ Sfc + x ~ Sfc ) G 

** n^fc(Sfc-^) 

Thus, without any restriction, we may assume that so, • • • , st are positive 
integers. We can rewrite the first equation in ( |2.17| ) as x 2 = 2wx — 1. From 
this relation, one can easily see by induction that 

x k = C/ fc -iHx - C4_ 2 H, for k = 1, 2, ... , (3.7) 

with the understanding that U-i(w) = 0. Using ([T?]), we can write the 
function in ( |3.1| ) in the form 

T 

fc =o Sfc HjVfcl s fc s jJ 
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where Q So ,...,s T (w) is the polynomial defined by Clearly, the sum in 
( |3.8[) belongs to Ay, so it remains to show that 

Qs ,...,s T {w)x e A v . (3.9) 



Using ( 2.17 ) we can write v in terms of w and x as 
v = (w - l) R (w + l) s (x - w), 

which gives 

(w - l) R (w + l) s x = v + w(w - l) R (w + l) s £ A v . 
Therefore, if T > max(i?, S), we have 

(w 2 - l) T x E A v . (3.10) 



The proof now follows from Lemma 3.3. □ 



4. Some identities satisfied by the Bessel functions 

The Bessel functions of imaginary argument are defined by the generating 
function 

Y^h{t)x k = e t{x+x ~ 1 V\ (4.1) 



whence Ik(t) = I_fc(t), and t) = (— l) fc // c (t). Differentiating (4.1) with 
respect to x and t, one gets 

kl k (t) = ~ (h-i(t) ~ h+i(t)) , (4.2a) 

and 

l' k {t) = \(h-i(t)+I k+1 {t)), (4.2b) 

respectively. Similarly, one can show that Ik(t) satisfy the modified Bessel 
equation 

(t 2 d? + td t -(t 2 + k 2 ))l k (t) = 0, (4.3) 

where dt = d/dt. The main result in this section is Proposition |4.l| below, 
which says that if q(j) is an odd polynomial of j, then the sum 

j>k 

j odd/even 

can be written as a finite linear combination of Bessel functions of the form 

finitely many j's 

where each coefficient (Xj(t) is a polynomial in t. Let us define a sequence 
of polynomials a" (t) for j G Z and n = 0, 1, . . . as follows 

«■$(,)_ S J0 ^ = {' /2 ifi = ° (4.4) 
2 I if j 7^ V ' 
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and 

a -+\t) = (t 2 d 2 + td t )a](t) + {fd t + |) (a] +1 (t) + a^t)) + 

t l(a] +2 (t)-2a](t) + a]„ 2 (t)). (4.5) 

From the symmetry of the defining relation, it follows that a"(t) = a" -(f). 
Moreover, it is clear that 

a™ = for j $ [-2n,2n], (4.6) 

and a"L 2n (t) = a% n (t) = t 2n+1 /2 2n+1 . For arbitrary j £ [-2ra,2n], e$(t) is a 
polynomial in t of degree at most 2n + 1 . 

Proposition 4.1. Ze£ fc and n be integers, n > 0. TTien 

2n+/c 

£ /«+i /j(t) = ^ a n_ fcWls(t)= £ a?_ fc (t)I.(t). (4-7) 

j>fc sSZ s=-2n+fc 

j=k+l (mod 2) 



Proof by induction. For n = we have to show that 

t 
2 



j=fc+l (mod 2) 



which easily follows from ( [4.2a| ). Assume that ( |4.7| ) holds for some n. Then, 
using the modified Bessel equation ([4.3|) we obtain 



1=0 



J2(k + 2l + l) 2n+3 I k+2l+1 (t) 

oo 

(f% + t $ - t 2 ) J^ik + 21 + l) 2n+1 I k+2l+1 (t) 

1=0 

(t 2 d 2 + td t -t 2 ) 52cZ_ k (t)i.(t), 



sez 



which shows that (4.7) holds for n + 1, upon using (4.2b). □ 

5. Heat kernel expansions on the integers 

In the case of the real line, the solution of the heat equation is not unique 
unless the class of solutions satisfies a condition of the form 

\u(x,t)\ < cie C2X ' 2 , ci,c 2 > 0. 

When one says that the Gaussian kernel is the fundamental solution of the 
heat equation one (implicitly) assumes that one is considering functions of 
moderate growth at infinity. Denote by u(n, m, t) the solution of 

Ut = L R:S u, (5.1) 
u\ t =o = 5 nm , (5.2) 



12 
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where L^s is the second-order difference operator (acting on functions of 
a discrete variable n £ Z) constructed in Section 2. We make the same 
implicit assumption here when we say that u(n, m, t) is the fundamental 
solution. 

From (2.15) and ( 2.21| ) and elementary spectral theory, it follows that 



u(n, m, t) 



r ( m ) e I 



r{m + 1) 2m Jq 



] p Tl (x)p m (x 1 ) — , (5.3) 

X 



where, for simplicity, we have omitted the dependence on the parameters 
r = (ri, r2, . . . ). Using ( 2.2C| ) we can write also the fundamental solution in 
terms of the wave and the adjoint wave functions as 



u(n,m,t) = ® e^ x+x ^ 

2m J c 



Pn{x)p* m+l {x)dx. 



(5.4) 



Theorem 5.1. The solution of (5JJ with initial condition (|5.2| ) can be 

written in the form 



u[n 



(5.5) 



finitely many j 's 



where (3j(n,m,t) are polynomials in t of degree at most 2T — 1, with T 
max(i?, S). 



Proof. From (|2~2C| ) it follows that 



r(m)r(n + 1) 

u(n,m,t) = — u(m,n,t). 

T(m + l)T{n) 

Thus, we can assume that k = n — m > 0. Around x = 0, we have the 
expansion 



Pn(x)p m (x 1 ) _ r{m + 1) 
x r(m) 



Pn{x)p* m+l {x) 



„n—m—l 



T(n)r(m) 



(r(n + l)r(m) +0(x)), 



which shows that 



res x=0 x j p n (x)p m (x 



-i^dx 



x 



for j > 1 — k. 



(5.6) 



For e = l,2 and i = 0, 1, . . . , T — 1 denote 

T 



q k+€+2i (j) 



J 



n 



3 



+ e + 2/) s 



k + e + 2* 11 (Jfe + e + 2i) 2 - (ib + e + 2/) 2 



(5.7) 



HEAT KERNEL EXPANSIONS ON THE INTEGERS 

<Ik+e+2i{j) 1S an °dd polynomial in j of degree IT — 1. We have 
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Ij(2t)x j +7 fc (2t)j 



j>l-k 



T-l 



+ E 

j>fc+2T 



E <lk+e+2i(j)x' 



(k+e+2i) 



8=0 



2 T-l 

+ EE 

e=l i=0 



E Qk+e+2i(j)Ij(2t) 



j>k+e+2i 
J=k+e (mod 2) 



-(fc+e+2i) 



(5.8a) 
(5.8b) 

(5.8c) 



where in the second sum, for every fixed j > k + 2T, we choose e = 1 or 2, 
so that j = k + e (mod 2). Denote by /j(x) the Laurent polynomial of x in 
the second sum in (fT^), i.e. 



T-l 

fj(x) = X~ ] - E <7fc+e+2i(j> 
i=0 



-(fc+e+2i) 



Notice that 



/;(*) = -J II^' 2 - ( fc + e + 2/) 2 )F S0iSl ,..., ST (x), 
z=l 

,(x) is the function defined by (|3.1[) with Sj = — (A: + e + 2i) 



where F S0)S1) ... )StV 

i = 0, 1, . . . , T — 1 and st = — j. Thus, by Proposition 3J, /j(x) S Ay, and 
therefore there exists a difference operator with support [—j, — {k + e)], 
satisfying 

fj(x)p n (x) = L f .p n (x). 
Since n — (A; + e) = m — e < m we see that 

(i.r 



fj(x)Pn(x)p m (x X ) = 



(5.9) 



From ( |5.3| ), ( ^ ) and ( |5.9[ ) it follows that the infinite sums in ( paj) and 
(5.8b) do not "contribute" to the fundamental solution u(n,m,t). Finally, 
the infinite sum in (5.8c) can be rewritten as a finite linear combination of 
Bessel functions with polynomial coefficients, according to Proposition 4.1, 
which completes the proof. □ 

We shall illustrate all steps of the proof by considering the case R = S = 1. 



Example 5.2. Let R = S = 1. From (2.2) and ( |2.3| ) one computes 
<j)i(n;r) = S-[(n;r) = n + n; 



ipi(n;r) 



n; r 



n 



+ J2(-2y- 1 jr J (-lfexp X> 2 )"'< 
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Denote for simplicity a = r\ and = Y^=2(~ ^9 3 r y The tau function is 
given by formula (|2,4|): 



r(n) 



n + a —n + a + (3 
1 2n + 1 - 2(a + 0) 



(5.10) 



The second-order difference operator L\ \ is given by formula 
d/dri = djda. 

From Q, (pl3|) and ( |2l8| ) we get 



5h, with 



x 



r(n)(x 2 — 1) 



n + a —n + a + (3 1 
n + 1 + a n + 1 — a — /3 x 
n + 2 + a -n - 2 + a + /? 



(5.11) 



From the last formula one can easily deduce that near x = we can 
expand p n (x)p m {x~ 1 ) as 



r(n + 1) 
r(n) 



(5.12) 



where 



7j 



(m — a — {3 + l)(n — a — + l)(n — m + j) 



r(n)T(TO 

+ (-l) J '(m + a + l)(n + a + l)(n - rn + j) 



(5.13) 



From ( p. 13 ) and ( |4.8| ) one can see that indeed u(n,m,t) is a finite linear 
combination of Bessel functions. Below, we shall illustrate how this can be 
seen following the proof of Theorem 5J. using just the first few coefficients 
in (|l2D. 

Proposition 3.1 tells us that Vs, I ^ 0, such that s = I (mod 2) 



x s x l 

6 Ay. 

s I 
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Following (|5.8|), we can write e t ( x+x ^ as 
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t(x+x~ x ) 



E Ij(2t)x j + I k {2t)x' 



j>l-k 



+ E ^ 

j>k+2 
j=k+l (mod 2) 

+ E ^ 2t ){ 

j>k+2 
j=k (mod 2) 



X-* - ^x-^ 



k + 1 



x i - ^—x-^ 



k + 2 



+ 



k + 1 



£ 

3>k 
j=k+l (mod 2) 



jU2t) 



x 



-(k+l) 



+ 



k + 2 



E iWt) 



j>k+l 
j=k (mod 2) 



X 



-(k+2) 



(5.14a) 
(5.14b) 

(5.14c) 
(5.14d) 
(5.14e) 



where k = n — m. Using (<L8), we can write the sums in (|5Tl4dD and (&U$) 



as 



E jlj(2t) = tl k (2t) tmd E jlj(2t)=tl k+1 (2t). 



j>k 
j=k+l (mod 2) 



j>fc+l 
j=fc (mod 2) 



Thus ( p.14 ) can be rewritten as 



j>i-k 

+ ^ 

j>k+2 
j=k+l (mod 2) 



x -j _ x -(k+l) 



+ £ 1,(2*) 



fc+1 



J_ x -(fc+2) 



j>k+2 
j=k (mod 2) 



k + 2 



+ x 



4(2*) + ^L^I^x' 1 + ] -^I k+1 (2t)x- 2 



(5.15a) 
(5.15b) 

(5.15c) 

(5.15d) 
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The sums in (5.15a), ( |5.15b| ) an d (|5.15c ) "do not contribute" to the integral 
( |5.3[) (see the proof of Theorem |5.1|). Thus 

%{2t) + r ^ T / fc (2t)x- 1 + j-^I k+1 (2t)x- 2 ^j 



u(n, m, t) 



e 2t TQS x= Q 



k + l 

T{m) Pn{x)p m {x^ 1 )- 



r(m + 1) 



x 



k+l 



r(m) /r(n + l) , , 71 . . 70 , . , 



T{m + 1) V r(n) 



fc + 1 



(5.16) 

where 71 and 72 are the coefficients in the expansion ( fOp . Using ( gig) 
we get the following explicit formula for u(n,m,t) 
„-2t 



n(n, m, t) 



r(m)r(n + l)/ n _ m (2t) 



r(m + l)r(n) 
+ 4i(/3 + 2a) (n + m - /? + 2)/ n _ m (2t) 

- 4t(2mn - f3n + 2n - /3m + 2m + /3 2 + 2a/3 - 2j3 + 2a 2 + 2)/ n _ m+ i (2t) 

(5.17)' 

If we put 5 = a + 1/2 and let /? — > 00 we get 



-2t 



n(n, m, i) 



T~mTn+lIn— m(2t) tl n — m (2t^ tl, t 



n— 771+I 



(2t) 



where r n = n + <5. Notice that this is exactly formula ( |l.6[) for the funda- 
mental solution computed in the introduction {R = 1, S = 0). 

The referee has raised the interesting possibility of using the formula 

/■t 



JAB 



l + A e sBA ds)B 







to get an alternative proof of Theorem |5.1| . This remains a challenging 
problem. 

Acknowledgments. We thank Henry P. McKean and the referee for 
suggestions that led to an improved version of this paper. 
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